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Using an approach based on quasiclassical Green's functions we present a theoretical study of 
transport in mesoscopic S/N structures in the diffusive limit. The subgap conductance in S/N 
structures with barriers (zero bias and finite bias anomalies) are discused. We also analyse the 
temperature dependence of the conductance variation 5S(T) for a Andreev interferometer. We show 
that besides the well know low temperature maximum a second maximum near T c may appear. We 
present the results of studies on the Josephson effect in 4 terminal S/N/S contacts and on the 
possible sign reversal of the Josephson critical current. 

1. Introduction 

The theory describing non-equilibrium transport in superconductors was developed 20 years ago (see, for example 
jij) and has been used to explain phenomena such as viscous flux flow, Josephson effects in superconducting weak links 
and the passage of the current over the interface between a superconductor and a normal metal (an S/N interface). 
Technological advances achieved in the last decade have enabled the fabrication of mesoscopic structures where 
the dimensions are less than the energy relaxation length and with well defined physical properties. In mesoscopic 
systems phase coherence is maintained and the distribution of quasiparticles may differ significantlty from equilibrium. 
Experiments on these structures have revealed a number of new effects such as, subgap conductance in SIN junctions 
(where I is an insulating layer), the oscillatory dependence of the conductance in S/N structures containing normal 
metal or superconducting loops, the non-monotonic dependence of the N film conductance on the temperature and 
voltage in S/N structures; long-range, phase-coherent effects, change of sign of the Josephson critical current in 4- 
terminal S/N/S structures when an additional current is driven through the N film. Although it is only recently 
that most of these effects have been explained using a non-equilibrium theory developed at the end of the 1970's (see 
||D). One method for studying these effects (ballistic systems in particular) is based on the Bogolyubov-de Gennes 
equations and the Landauer formula for the conductances In this paper we use another approach based on a 

microscopic, quasiclassical Green's function technique, which has been successfully applied to systems with a small 
mean free path (i.e. diffusive regime), to analyse these effects. The method of matrix, quasiclassical Green's functions 
is a convenient and powerful method for studying transport in mesoscopic S/N structures Q. The quasiclassical 
approximation means that all the Green's functions are spatially averaged over distances of order of the Fermi wave 
length Pp 1 . This is equivalent to the intergration in momentum space over the variable £ p = (p — pp)vp. To describe 
transport and non-equilibrium properties we need to know the retarded (advanced) Green's functions G R ^ and the 
Keldysh function G. In the case of superconducting systems each of these functions is a 2x2 matrix elements of which 
describe normal excitations and a condensate. These matrices are grouped together, defining the 4x4 matrix Green's 
function G as 



where G R ^ are the matrix retarded (advanced) Green's functions and G is the Keldysh function. The functions 
G R ( A ) describe thermodynamical properties of the system such as the density of states (DOS), condensate currents 
and so on, while the function G describes the kinetic and nonequilbrium properties of the system. The matrix G is 
related to a matrix distribution function / by 

G = G R f - fG A (2) 

Where the elements of the matrix / = f l + fa z describe particle and hole-like excitations. The component f Q is 
used in calculating the order parameter and supercurrent, while / describes the branch imbalance and determines the 
spatial distribution of the electric field (see for example ) . The matrices G R ^ can be presented in the form, 

qR(A) = G R(A) dz + pR(A) (3) 
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In a bulk superconductor G R ^ = e/^ (A) , F R ^ = F R ( A 1 (ia x sin X + iv y cos X ) , F R ^ = A/^ (A) , ^ {A) = 
y/ (e ± ir) 2 — A 2 , r is the damping rate in the excitation spectrum of the superconductor and \ is the phase of 
the order parameter. The current in the system is expressed in terms of the function G 

I = (crd/8)Tra z J de(G R d x G + Gd x G A ) (4) 

This expression is written for the one-dimentional case, where d is the thickness of the system. In the diffusive limit 
the supermatrix G in the N film of a S/N mesoscopic system obeys the equation 

Dd x [Gd x G] +ie[a z ,G] = 0. (5) 

where D is the diffusion coefficient and a z is the Pauli supermatrix. We neglect the inelastic collision integral by 
assuming that the length of the N film L is short enoughierh >> [^ lT e-eT^ lT e-ph\^ { e Th = h,D / L 2 is the Thouless 
energy and r e _ e ,r e _p^ are the electron-electron and electron-phonon collision times). To solve Eq.(|^) we must have 
the complementary boundary condition, 

D(Gd z G) =(e b d N ) [G,G S ] (6) 

the z-axis is normal to the plane of the S/N interface. The boundary conditions for the quasiclassical Green's functions 
G were derived for the general case by Zaitsev Q and reduced to the simple form (^) by Kupriyanov and Lukichev 
in the dirty case. The applicability of Eq.(||) is discussed in Ref. [||. The distribution functions in the reservoirs 
fo and / are assumed to be in equilibrium and for example, for the structure shown in Fig. 1(b) they have the form, 

f (±L) = F oN (e) = [tanh((e + eV)0) + tanh((e - eV)0)}/2, (7) 
f(±L) = F N (e) = [tanh((e + eV)/3) - tanh((e - eV)/3)}/2, (8) 

where (3 = (2T) _1 and ±V are the electric potentials in the reservoirs at x = ±L (we have assumed the system is 
symmetrical). In order to find the distribution functions f(e,x) in the N film, we consider element (12) of Eq.(|^) 
(i.e. we consider the equation for the Keldysh matrix G). Consider for simplicity a strucutre of the type in Fig. 1(b) 
with one superconducting strip in the centre of the N film. Multiplying this equation by a z and taking the trace, we 
obtain, 

d x {M{x)d x f) = (9) 

where M(x) = (1/2) (1 — G R G A — F R F A ) and we have assumed that the width of the S strip (w) is small compared 
to L. In this limit the distribution functions /(e,0) and /s(e, 0) are equal to zero (V(0) = 0). Integrating Eq.(g) and 
taking into account boundary condition (@), we obtain, 

J={F N /L)<M- 1 >, (10) 
here < M _1 >= dxM^ 1 (see ||) and J(e) is the integration constant which determines the current, 

/>OC 

I = (ad/2) / cfeJ(e). (11) 
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Therefore, the problem of finding the I(V) characteristics is reduced to finding the retarded (advanced) Green's 
functions which determine the function M(e, x). Using the well known substitution, G^" 4 ' = cosh and F R ( A > = 
sinhu R ( A \ the equation for the retarded (advanced) Green's functions G R ^ may be presented in the form, 

d xx u R{A) - (k R{A) ^ smhu R ^ A) = -k 2 w(5 (x + L x ) + S (x - L x )) [f* {A) cosh u - G R S {A) sinh «] (12) 

We have used Eqs.(||) and @ and assumed that w « L, where fc 2 = e&/ D,(k R ^) 2 = 7 =p 2ie. The term on the 
right hand side of Eq.(|l2|) is the source of the condensate (the proximity effect). Solutions to Eq.(|l2[) can be written 
in an explicit form in some limiting cases. For example in the case of a weak proximity effect << 1) Eq.(|l2]) 

can be linearised (i.e. sinhu^ w u fl ), then a solution of Eq.(p^|) is easily found. This case occurs if the ratio of the 
resistance of the N film i?jy and the S/N interface resistance Rs/n 1S small, r = Rn/Rs/n << !• The comparison 
of an exact numerical solution and the solution of the linearised Eq.(|l^) shows that even when r « 1 the difference 
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between these solutions is less than 10%. Another limiting case corresponds to a short N film (exh >> (T,V)), now 
Eq.(|l2|) can be averaged over the length jlO 11 1. We use this procedure to analyse the subgap conductance of the 



structure shown in Fig. 1(a) 

2. Subgap conductance of S/N/N' structures 

Considering the S/N/N' system shown in Fig. 1(a) and assuming the middle N film is short, i.e. cth = D/L 2 >> 
[T, eV] ( L is the length of the N film and V the applied voltage), we can average Eq.(|l2|) over the length using the 
boundary condition (||) and obtain for F R ( A ^ and G R ^ , 

f r(a) = t9 ^ ■ G R ^ = e±%1 2 (13) 

Y (e ± ilf -£g y (e ± «7) 2 - £g 

Here e 9 = pD / (2RbsnL), 7 = pD/(2RbNnL) and we have assumed that A >> (e,e g ). Rbs.Nn are the interface 
resistance per unit area. The energy e g is a minigap in the excitation spectrum of the N film induced by the proximity 
effect [ |12|Jl3]| . The damping constant 7 is due to the N/N' contact, which increases with decreasing the N/N' interface 
resistance. In Fig. 2(a) we show the energy dependence of the DOS in the normal film for different 7, 

u N = [G R - G A )/2 (14) 

We see that if 7 >> e g the DOS in the N film for T < T c is close to the DOS when T > T c , and a subgap arises 
when 7 << e g . 

In order to calculate the conductance of the structure, we need to find the distribution function / from Eq.(^|), 
taking into account boundary condition ([3]). Carrying out some simple calculations, we obtain the dimensionlcss 
conductance ||, 

f°° 1 n-r^r rir- 1 , 1 / ^-Th \ / ^g \ , ^Th 



AT C°° 
S = RN dV = J d€ P F Nl M 



(— + (15) 

v d £ g 7 Ae g 



Here F' N = dF N /d(eV), M is defined in Eq.@ and 

1 



A = -[(G R G )(Gg - Gg) - (F R + F )(F R + F#)] (16) 

The terms in the square brackets in Eq. fllq) determine the normalised resistances of the N film, N/N' and N/S 
interfaces respectively. The first term in Eq.(|16[) stems from the usual quasiparticle current, which is non-zero when 



eV > A. The second term in Eq.(16) represents the subgap conductance due to Andreev reflection processes and is 
non zero when e < A (if the damping T in the superconductors is negligible). This term also leads to the so called 
interference current in Josephson tunnel junctions. In our case this term arises due to the proximity effect inducing 
a condensate in the normal film. In Fig. 2(b) we plot the dependence of the zero-temperature conductance S verses 
normalised voltage, S — G n (dl /dV), where G n is the normal state conductance. The parameter eTh/tg, the ratio of 
the S /N interface resistance to the N film resistance is chosen large so that the resistance of the system is determined 
by the interface resistance. If the damping rate 7 is large in comparison with e g , then the DOS in the normal film 
is only weakly disturbed by the proximity effect and the peak arising in the subgap conductance is located at zero 
voltage. The S/N interface resistance dominates in this case, and the subgap conductance is caused by the interference 
current (the second term in Eq.(|l6|)), or by Andreev reflection processes. In the opposite case of a small damping 
(7 << e g ) or a large N/N' interface resistance a peak in the DOS of the N film appears in the voltage dependence 
of the subgap conductance (curves for 7 = 0.1 and 7 = 0.5). The peak in S(V) is caused by a peculiarity in the 
quasiparticle current at eV = e g (the second term in Eq.(fl6|)). 

It is interesting to note that the minigap in the N film e g may be modulated by an external magnetic field if 
two superconductors (instead of one) are used in the structure (see for example (m) ) , then e g depends on the phase 
difference between the superconductors <fr, e g = e g o\ cos(f>\. We see from Fig. 2(b) that increasing 7 with respect to e g 
may lead to an increase in the conductance for some voltages. This means that an applied magnetic field H ((f> ~ H) 
may cause an increase in the conductance. This effect was predicted for a two-barrier structure in refs 0,^5) and 
observed experimentally ^5),^6|. Note that the subgap conductance was first observed in a SISm system (Sm is a 
heavily doped semiconductor) by Kastalsky and Kleinsasser |3(| and later by others |3^-Q , and theoretically studied 
in many works P,[36"|-j39[ . 

3. Nonmonotonic temperature and voltage dependence of the conductance. 

In this section we analyse the voltage and temperature dependence of the conductance for the structure shown 
in Fig. 1(b) (Andreev interferometer). Since the late 1970's it has been known that the conductance (G) of S/N 
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mesoscopic structures depends on temperature (T) (and voltage (V)) in a non-monotonic way (see reviews 
This behaviour was first predicted in Ref. O] where a simple point S/N contact was analysed. The authors of Ref. 
p6| , using a microscopic theory and assuming that the energy gap in the superconductor (A) is much less than 
the Thouless energy exh = fiD/L 2 , showed that the zero-bias conductance G coincides at zero temperature with its 
normal state value (G n ). With increasing T, G exhibits a non-monotonic behaviour, increasing to a maximum of 
Gmax ~ 1-25G„ at T m rj A(T m ) and then decreasing to G n for T > T m . 

Recently mesoscopic S/N structures have been fabricated in which the limit A >> em is realised. In this case 
Nazarov and Stoof JlTj (also see flg|-|20[| ) argued that the temperature dependence of the conductance G has a 
similar non-monotonic behaviour with a maximum at a temperature comparable with the Thouless energy, while 
simultaniously Volkov, Allsopp and Lambert |21| predicted that the voltage dependence of the conductance in an S/N 
mesoscopic structure (Andreev interferometer) has a similar form with a maximum at eV m ~ £tk- This non-monotonic 
behaviour has been observed both in very short S/N contacts [^2j and in longer mesoscopic S/N structures p3|-p5| , p7|] . 

Considering the system shown in Fig. 1(b) and assuming a weak proximity effect (IF^^I << 1), Eq.(|12[) and 
< A/ -1 > can be presented in the form, < M -1 >= 1— < to > ^0|, where 

(m) = -Trl(P R ) 2 + (P A ) 2 - 2F R F A \ (17) 



The first two terms in Eq.(|17]) determine a change in the DOS of the N film due to the proximity effect. This can 
be seen from Eq.([u]) by taking into account the smallness of F R ^ and the normalisation condition Q 

(G R Wfi + (F R W) 2 = l. (18) 

The DOS in the N film is reduced at small energies by the proximity effect and therefore these terms reduce the 



conductance of the system. The last (anomalous) term in Eq.(17) is analogous to the so called Maki-Thompson term 
in the paraconductivity and increases the conductance (see Fig. 3). At zero energy both contributions cancel each 
other. This is seen from Eq.(|l2]). At e = the equations for u^- 4 ) coincide with each other and F R — F A , hence 
< to >= (see Eq.([l7|)). In order to calculate the variation of the normalised conductance SS, one has to solve 
Eq.(|l2[) and to find the functions F R ( A \x). In the case of a weak proximity effect Eq.(|l2|) can be linearised, and a 
solution is readily found, 

F R(A \e,x) = iv y rF R[A) \smh{k R{A) {L - x)j\ cos(0/2)/((9cosh6O 2 (19) 

where r is the ratio of the N film resistance and the S/N interface resistance in the normal state, is the phase 
difference between the superconductors, for simplicity we assume that L\ « L. Knowing F R ( A > we can easily 
calculate < to > and SS using Eqs.( ^0[]Tl| ). The function < to > which determines the voltage dependence of the 
zero-temperature conductance is equal to, 



< m >- 



Re /sinh^ _ \ /{dcQshe)2 /sinfi^ _ ( 



\ 28 J V 2 °i 26) 2 



(l + cos0) (20) 



where = 6\ + 162- The first term in Eq.([20j) is related to a change in the DOS and the second term is due to the 
anomalous term < F R F A >. In Fig. 3 we present the total conductance variation SS vs the normalised voltage V 
(in units of the Thouless energy) at zero temperature, the contribution SSdos due to a change of the DOS, and the 
anomalous term SS an for = 0. The variation SS has a maximum at eV « cth |2l}| . The temperature dependence 
(^(T) found with the help of Eqs. ( |iof]ri| ) has a similar form provided T « A, er/, << A p7|-p0[. It is interesting 
to note that as T increases 8S an decays slowly (non-exponcntially) as (eTh/T) at T » tTh- This law leads to 



long-range, phase-coherent effects in mesoscopic structures p0| , p8 -30|. One can see from Eq.(|19|) that the functions 
\pR{ A ) J have a peak at e w A because the functions F R ^ have a singularity at the energy |e| = A (see Eq.(|J)). With 
increasing temperature T the order parameter A(T) starts decreasing and near T c may be comparable with 67^ , hence 
the function < to > and SS may have one more maximum. In Fig. 4 the temperature dependence of SS is presented 
over a wide temperature range (0 < T < T c ) for a structure similar to that shown in Fig. lb (where L\ = 0, i.e. only 
one superconductor is in contact with the normal film). We see that besides the main, low-temperature maximum in 
SS (T) there is an additional maxinum near T c . The magnitude of both maxima depends strongly on the depairing 
rate 7 in the N film (in order to account for 7, the energy e must be replaced with e ± ij) [j35f . 

4. Negative Josephson current in a 4 terminal S/N/S structure. 

In the preceding section we calculated the conductance S of the Andreev interferometer (see Fig. lb) , determined 
by the nonequilibrium distribution function / (see Eq.(^])) which varies in space. In the limit of a weak proximity 
effect we obtain from Eq.(||) 
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f(e,x) = J(e)[x+ / dx 1 m(x 1 )] 
Jo 



(21) 



Here J(e) is a function even in energy e and defined in Eq.(lO). Therefore f(e,x) increases almost linearly from 
zero (at x = 0) to F/v (at x = L). The function < m > determined by the condensate functions depends on the phase 
difference <p (see Eq.(20)), increasing (for example, in an applied magnetic field) causes the function < m > and 
the conductance S to oscillate. This oscillatory behaviour of the conductance has been observed in many experiments 
IIII^JHHHQ startin g from tne pioneering work Q and later explained theoretically . 

One can formulate the inverse problem, how the non-equilibrium distribution function affects, for example, the 
condensate current, which is determined by the expression, 



h = —Tr 



de F R dF H - F A dF 



So 



(22) 



where f is a function odd in energy e and determined in Eq.(^). By taking the trace of this equation in the limit 
of a weak proximity effect we obtain, 



d x [{l - m )d x f ] = 



(23) 



here m Q = (l/8)Tr(F R + F A ) 2 . The current I s is of the order of the small parameter r 2 , therefore in finding f D 
we can neglect m and obtain that f Q is constant in space and equal to / (e) = -Fjvo (see Eq. (?])). The function 
Fno{c) is an equilibrium distribution function (odd in e) in the N reservoirs at a given potential V, and is shifted 
with respect to the function Fqs = tanhe/? by the value V (the potential in the superconductors is set equal to zero). 
The function Fnq is an even function of V and therefore does not depend on whether the additional (control) current 
flows from the N reservoirs to the superconductors (the potential V has the same sign in both reservoirs) or from one 
N reservoir to another one (V has opposite signs in the reservoirs) . A relative shift of the distribution functions in N 
and S controlled by the applied voltage V leads to a dependence of the Josephson critical current I c on V. In order to 
determine this dependence, we need to find the condensate functions F R ( A \ obeying the linearised Usadel equation 
which in this case has the form, 



[S (x + Li) + 5{x- Lx) 



(24) 



the functions Fg are given by Eq.(||). Solving Eq.(|24]) and substituting the solutions into Eq.(p^), we find 
/, = I r sin (A, where 



I C {V)R L = (7rT)r 2 3ff 



E 



sinh 2 (0 n (£-£!)/£) 



— i A 2 + (w n + ieV) 2 



9 n sinh(26»„) 



(25) 



here 9 n = L[(2(uj n + ieV)+j)/D] 1 / 2 , uj n — nT(2n + l) is the Matsubara frequency. In Fig. 5 we show the dependence 
of the Josephson critical current I c on V. We see that I c changes sign when eV w exh- The I c sign reversal effect 
(it - contact) was first predicted in Ref (i^] where a S/M/S contact was considered (M is an insulating layer with 
magnetic impurites). The same effect was analysed later in Ref Q] where a S/F/S contact was considered. The 
possibility to change of the I c sign by a control voltage (or current) in a 4-terminal S/N/S diffusive structure was 
first studied in Ref |45| for the 1- dimentional case of a short spacing between superconductors (L\ tdis) and in Refs. 
pp| , p9 46 4?|] for the case of a longer spacing. The physical reason for the sign reversal of I c in 4 terminal S/N/S 
structures is a difference in the distribution functions f in the N film and in the superconductors. These functions 
have an equilibrium form corresponding to the electric potential V and zero respectively. The 1-dimensional case 
was analysed in [^o| , ^9||4^j46]] . In Ref jl?]] the 2-d case was analysed, however in finding the distribution funtion f D 
and the condensate functions F R ^ A ^ the authors considreed opposite limits, d s « d and d s » d, where d s ,d are 
the width of the S and N films respectively. Therefore the region of applicability of the results obtained in Ref |47]] 
remains unclear. One can show that qualitatively Eq.(|25|) in limit d s « d remains valid for the 2-d case, i.e. I c Rn 
is proportinal to the small parameter r 2 (r is the ratio of the N film resistance to the resistance of the S/N contact). 
Similar effects in ballistic systems were studied in Refs. 49 - 5l|] . 

Note, a formal analogy between the two possibilities of the I c sign reversal effect: a structure with magnetic 
material and a shift of the distribution functions in a 4-terminal S/N/S structure. In the first case the Green's 
functions F R ^ (e) in the F layer differ from those in the N layer by a shift in energy e — > e + h, where h is the 
exchange energy 44 1. The distribution function f a in the F layer coincides with the distribution function in the 
superconductor Fqs = tanhe/3. Meanwhile in the 4-terminal S/N/S contact the distribution function in the N layer is 
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shifted with respect to Fso, /o(e) = [Fos(e + eV) + i*bs( e — eV)]/2. Making the substitution e — > e ± eV in the integral 
in Eq. (|25|) , one can reduce the case of a 4-terminal S/N/S contact to the case of a S/F/S contact with h = eV. The 
sign reverse of I c in S/F/S contacts has, presumably, not yet been observed. On the other hand in a recent publication 
|tt8| experimental evidence for the change of the sign of I c in 4-terminal S/N/S contacts was obtained. 

5. Conclusions. 

The quasiclassical Green's function technique we have used complements the approach based on the scattering 
matrix method and the Landauer formula for the conductance. This technique allows one to describe the transport 
properites of S/N mesoscopic structures over a wide range of parameters. With the help of this approach the conduc- 
tance of a double-barrier S/N/N' structure was calculated. We have shown that in the differential subgap conductance 
S has a peak situated at zero or finite bias voltage (depending on the relation between the S/N and N/N' interface 
resistances). When the temperature is lowered below T c , the variation of the conductance may be both positive and 
negative. Therefore the conductance decreases or increases with increasing an external magnetic field. 

The conductance in S/N structures of the types of the Andreev interferometer or S/N point contact depends on T 
or V in a non-monotonic way. As a function of temperature T the zero-bias conductance variation 6S increases from 
zero, has a maximum at T% of order the Thouless energy exh (if (-Th << A) and decays to zero at T — > T c . Near T c a 
second smaller peak in 8S may occur at a temperature T2 corresponding to the condition A(T2) ~ £th- 

An interesting effect may occur in 4 - terminal S/N/S mesoscopic structures (see for example Fig. lb). If an additional 
(control) current flows through the N film, then the distribution functions in the N and S films are different. Due 
to this difference the critical Josephson current may change sign and the S/N/S contact may become a 7r-contact. 
Most of the effects discussed in this paper have been observed experimentally, although some predictions such as the 
appearance of the second maximum near T c in 5S have not yet been studied experimentally. 

A. F. Volkov is grateful to the Royal Society, to the Russian grant on superconductivity (Project 96053) and to 
CRDF (project RP1-165) for financial support. 
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FIG. 1. The structures considered. 



(a) 



(b) 





FIG. 2. (a) The energy dependence of the DOS for the structure in Fig. la, for different 7; here A = 20, T = 0.1 and e<j = 10. 
(b) The zero temperature conductance S (Eq.(|l5|)), normalised to the conductance of the normal metal in the norml state, 
dependence on voltage, using the same parameters as for Fig. 2a. All energies are measured in units of erh- 
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FIG. 3. Normalized conductance SS vs. normalized voltage eV/e^ at T = and vs. normalized temperature T/sl at V = 
for the structure shown in Fig. lb (Li = 0). 




0.5 

T/Tc 

FIG. 4. Dependence of SS on temperature, for the cross geomentry in the weak proximity limit, for different depairing rates 
7. With Ao = 20, T = 0.1,Li/L = 0. The inset is an enlargement of the second maximum around T/T c = 1.0. Note the 
5S = SS{2r b ) 2 . 




FIG. 5. Critical current 7 C dependence on voltage for the 4 terminal Andreev interferometer, for several temperatures (T). 
The critical current is measured in units tTh /(eR L r 2 ) and T, V in £Th;A = 20, 7 = 0.2, Li/L = 0.5. 
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